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Abstract. We compute the d-dimensional critical exponents corresponding to the wave function 
and mass renormalization of the quark in QCD in the Landau gauge at a new order, 0(1/Nj), in 
• the large Nf expansion. The computations are simplified by the establishment in d-dimensions 

of the critical point equivalence of QCD and the non-abelian Thirring model beyond leading 
order. The form of the 0{1/Nj) coefficients in the MS quark mass anomalous dimension at five 
$_i ■ loops is deduced and compared with the numerical asymptotic Pade approximant prediction. 
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The renormalization group equation, (RGE), plays a fundamental role in our understanding 
of the properties of quantum field theories. Central to this equation are the renormalization 
group functions such as the /^-function and the anomalous dimension, 7(5), and for theories 
involving massive fields the mass anomalous dimension, j m (g), also appears. Ordinarily such 
functions are computed order by order in the loop expansion. Though this clearly becomes 
more difficult at successive orders due to the increase in the number of Feynman integrals which 
need to be computed. Despite this problem it has been possible to determine the /3-function 
and quark mass anomalous dimensions in quantum chromo dynamics, (QCD), at four loops in 
the MS renormalization scheme, [1-11]. Such calculations represent a significant achievement 
especially given that the order of 10 4 Feynman diagrams have been evaluated at four loops to 
deduce the renormalization constants. To proceed further in perturbation theory will clearly 
be a colossal if not impossible undertaking. However, there do exist methods which probe the 
perturbative structure of the RG functions from another direction. Given that they depend not 
only on the coupling constant, g, but other parameters in the theory such as the dimension of 
any internal symmetry groups, one can equally expand in one of these other variables. This 
will still involve the computation of Feynman integrals but not those associated with ordinary 
perturbation theory. In QCD, with Nf flavours of quarks, they can be expanded in the large 
Nf expansion where 1/Nt behaves as a bona fide perturbation parameter. Explicit details of 
the technique have been recorded in the literature but we emphasise that the extension of the 



O(N) a model methods of Vasil'ev et al, [jC], 13|, to four dimensional gauge theories, [14- 
18], offer the best and most efficient strategy to computing 1/Nf information in QCD. For 
example, the anomalous dimension of the twist-2 operator dimensions which are fundamental 
to the operator product expansion used in deep inelastic scattering have been computed at 
0(1/Nf) and to all orders in the strong coupling constant, [jl9(| . These results have been crucial 
in confirming the correctness of the explicit perturbative 3-loop results of [^] which play a 
key role in the full two loop evolution of the QCD structure functions. Clearly, given this 
important overlap with the current activity in explicit perturbative calculations, it is crucial 
that the large Nf method is developed to the next order, 0(1/Nf). However, before such a deep 
inelastic programme can proceed, various fundamental computations need to be performed. In 
any ordinary renormalization the wave function renormalization constant is always computed 
first before determining the renormalization constants of the other parameters of the field theory. 
Likewise in the large Nf programme, the wave function critical exponent, 77, must be deduced 
prior to the calculation of any other operator or field dimension. Through the critical RGE, 77 is 
related to j(g) by 77 = 2j(g*) where g* is the value of the critical coupling at the c?-dimensional 
fixed point of the QCD /3-function. 

Whilst 77 has been determined at 0(1 /Nj) in an arbitrary covariant gauge as a function 
of d in QED, [16|, the extension of that calculation to QCD has not yet been provided. This 
is one of the main aims of this letter where we will determine 77 in a simpler model in the 
same universality class as QCD, known as the non-abelian Thirring model. As we will argue 
it has a simpler form than QCD since the triple and quartic gluon vertices which are present 
in Yang-Mills theories are absent but the expression we obtain for 7/ will correspond to the 
quark wave function renormalization constant in QCD itself. Clearly this reduces the number 
of Feynman diagrams which need to be considered. Indeed the 0(1/ Nf) calculations of say, [JOJ, 
were performed in the non-abelian Thirring model, (NATM). The connection with this model 
had previously been investigated at leading order in the large Nf expansion in [21|. Although 
the provision of 77 is fundamental to any future 0(1/Nf) calculation, we have also computed the 
critical exponent which relates to the quark mass dimension, 7 m (g*), at 0(1/Nf ). There are 
several reasons for carrying out such a calculation. First, 7/ is gauge dependent and is therefore 
not a fully meaningful physical quantity. On the contrary r ) m (g*) is known to be a gauge 
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independent (and scheme independent) exponent. This provides us with a non-trivial check on 
our computation, aside of course from the comparison with the explicit four loop MS results of 
JToj , |IH |. In addition, the Feynman integrals with a [tpijj] insertion are closely related to those 
required to deduce r/ and so their computation does not in fact represent a significant amount of 
extra work. Moreover, as the calculation of r y m (g*) involves the insertion of the mass operator 
[ij)ip\ m a 2-point quark Green function, it very much lays the foundation for the computation 
of the dimensions of similar operators which are bilinear in the quark fields such as those which 
arise in deep inelastic scattering. 

The QCD Lagrangian in d = 4 — 2e dimensional Euclidean space reads 
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where is the quark field belonging to the fundamental representation of the colour group, 
1 < I < Nf, is the gluon field, c a and c a are the ghost fields in the adjoint representation of 
the colour group, £ is the covariant gauge parameter and g is the coupling constant. The field 
strength tensor i* 1 " and the covariant derivative are defined as F® v = d^A a v - d^A* + f abc A a ^A c u 
and = (d^ — iA^T a ), where T a are the group generators in the corresponding representation 
and f abc are the structure constants with [T a ,T b ] = if abc T c . To ensure the coupling constant, 
g, is dimensionless below four dimensions we rescale it in the standard way by setting g — > M e g, 
where the parameter M has dimensions of mass. 

The partition functions of the theory defined by the Lagrangian (P read 



(Oi(x! 



...On(xn)) = Z~ l J D<S> 1 (x 1 )...O n (x n )exp{-5} , 



where <3? = {A, ijj, ip, c, c} is the set of fundamental fields and Oi{xi) represent a basic field or 
a composite operator. The divergences arising in the calculation of (||) are removed at each 
order of perturbation theory by the renormalization of the fields and parameters entering the 
QCD Lagrangian, as well as by renormalization of the composite operators. Hence the renor- 
malized 1-particle irreducible n-point Green functions with the insertion of k multiplicatively 
renormalizable composite operators satisfy the RGE 



Md M + Pgdg + - n$7$ + ^70 4 r(xi, . . . , x n+k , M, g,£) = , 
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where 70,(5) is the anomalous dimension of the operator 0,, n$7$ = n^7A + nf^ip + "-c7c and 
the beta functions of the couplings g and £ are given by (3 g = Mdug and /3g = M8mC- It should 
be noted here that the correlation functions of gauge invariant objects do not depend on the 
gauge fixing parameter, £, and in this case the term fi^d^ drops out of Eq. (||). 

Our analysis relies heavily on the existence of a non-trivial infra-red, (IR), stable fixed point 
g* of the d-dimensional /3-function, P g (g*) = 0, for large values of Nf. The /^-function has been 
calculated in MS using dimensional regularization and in the notation of the [22] with a = 
(^^isiUll 



Pa (a) 



(d 



4)a + 

11 

72 



-T F N f 



11 

IT 



C A 



a 2 + 



^C F T F N f + ^C A T F N f 



17 
12 



C F T 2 N 2 + ^C A T 2 Nf + ^Cj,T, \ 



1415 ^ 2rrAT , 2857 ^ 



864 



■C A T F N f + 



1728 



a 



16 

4 + 0(a 5 ) , 



205 
288 



C f CaT f Nj 



(4) 



3 



from which it follows that 



a* = j^r + ^^(33C A e-[27C F + 45C A ]e 2 + 0(e 3 )) + O (±j . (5) 



The Casimirs for a general classical Lie group are defined by 

Tr (T a T b ) = T R 5 ab , T a T a = C F I , f acd f bcd = C A S ab . (6) 



It immediately follows from Eq. (||) that the Green functions of gauge invariant operators are 
scale invariant at the critical point g*. In other words G(Xxi) = X Di G(xi), where Di is the 
scaling dimension of the corresponding Green function. Moreover, due to the IR nature of the 
fixed point, this index determines the power of the leading term of the IR asymptotic behaviour 



of the Green functions (pi — > 0) near the critical points [23]. On the contrary Green functions of 
gauge dependent objects, such as the propagators of the basic fields which will in general depend 
on £, are not invariant under scale transformations. Although one may restrict attention from 
the outset to gauge independent quantities, since they have physical meaning, it is also possible 
and convenient to choose £ so that all Green functions are scale invariant. Evidently, this is 
equivalent to the condition /%(<?*,£*) = 0. Since 

Pt(j9,S) = -2Z(e + lA + P g /g) , (7) 

one concludes that the equation /%(<?*,£*) = has two solutions. One is = whilst the other 
is 7,4(5*; £*) = — £■ Bearing in mind that our main aim is the development of the 1/Nf expansion 
we choose the first solution, £ = 0, since the latter gives £ ~ Nf, which leads to problems in 
the construction of the 1/Nf scheme. The origin of the above two solutions for £ becomes more 
transparent if one tries to write down the most general form of the gluon propagator satisfying 
the requirements of both scale and gauge invariance. Indeed, scale invariance yields 

M 2e / \ 

g ^(p) = -^{ AP ^ + BP H > (8) 

where P^ v and pj\ u are the transverse and longitudinal projectors, respectively, and A and B 
are constants. As is well known, j^3|, radiative corrections do not contribute to the longitudinal 
part of gluon propagator. Hence, g\ v = ^g 2 M 2e pjl l ,p~ 2 . This implies, that if a 7^ 1 then £ must 
vanish, £ = 0. On the other hand for £ 7^ then one must have a = 1 which is easy to check 
is equivalent to 7^ = — e corresponding to the canonical dimension of the field. Earlier work 
concerning the relation of scaling and conformal symmetry in the context of gauge theories has 
been given in [24|. 



It is well known from the theory of the critical phenomena [23| that the critical properties 
of the system do not depend on the details of the interactions but is determined mainly by 
"global" characteristics such as symmetries and the dimension of spacetime. Thus different 
systems may exhibit the same behaviour at the critical point. An example of this universality is 
the fixed point relation between the Heisenberg ferromagnet and 4> A field theory. In what follows 
we construct the theory which belongs to the same universality class as QCD but which has a 
simpler structure. We first develop the 1/Nf expansion for calculating correlators of the type 
given in (||) . This can be achieved in the standard manner by integrating over the fermion fields 
in the functional integral which yields the following effective action for the gauge field, 

rr / M~ 2e M~ 2e \ 

Sf = Nf ( -trlnO - i^T a ) + ^r{F%? + ~^ri dA f I + V ( D " C T ( 9 ) 
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where bearing in mind that g 2 ~ 1/iVf we have set g 2 = g 2 /Nf. The evaluation of the functional 
integral with action (^) by the saddle point method generates the systematic expansion for the 
correlators. If we now examine the IR asymptotic behaviour of the Green functions in this 
approach and first of all consider the gluon propagator to first order in 1/Nt in the Landau 
gauge, £ = 0, it is 



d _ 2 M- 2 " 2 
^9 



, (io) 



where the first term in the brackets arises from the fermion loop. It is obvious that the contri- 
bution coming from the (F^ u ) 2 term is less singular in comparison with the contribution of the 
fermion loop in the asymptotic limit p — > 0. It can be easily checked that the diagrams with 
vertices contained in {F^ v ) 2 also do not contribute to the leading order of the IR limit. So one 
concludes that this term does not influence the critical properties of the theory and according 
to the general scheme should be excluded from action. Therefore we obtain the theory defined 
by the action 

S = iP{@-i4 a T a )iP + ^(U- e / 2 dA) 2 + d^cfd^c + f abc d»cAlc c , (11) 

which in the Landau gauge has the same critical behaviour as QCD. Of course, for gauge 
independent quantities it is true in any gauge. We have also modified the form of the gauge 
fixing condition, in order that the transverse and longitudinal parts of the gluon propagator have 
the same momentum dependence. In fact, to derive (|ll|) one can start from the theory with 
manifest gauge invariance which is determined by the action S = — iA^ a T a )^ with ghost 
and gauge fixing terms in turn arising from the application of the Faddeev-Popov procedure 



to the functional integral. Power counting shows that (11) is renormalizable within the 1/iVf 
expansion. Of course, we assume that a gauge invariant regularization is used. The renormalized 
action then takes the form 

S R = Z x $ ^ - iZ 2 ^4 a T a 4, + ^L(n-^ 2 dA) 2 + Zad^&c + Z^d^cA^c . (12) 

Due to the Slavnov- Taylor identities the renormalization constants Zj are related by 

ZsZf 1 = Z 4 Z 3 ^ 1 . (13) 

This was used in the exponent formulation to determine the ghost anomalous dimension at 
0(1/Nf) in p5| . It should be noted that in the Landau gauge Z4 = 1. As was proved above 
in the Landau gauge, the critical properties of QCD and this new theory which we shall refer 
to as the non-abelian Thirring model are identical. Therefore one can use the NATM model 
to deduce the QCD RG functions. This equivalence at leading order in 1/Nf was noted in pi] ] 
and used to deduce various exponents at 0(1/N f ), 0, H, 0, 0, HI, H- The extension of 



these calculations to 0(1/N 2 ) requires special care. The main one is the necessity of using a 
gauge invariant regularization which was not crucial at 0(1/Nf). The conventional dimensional 
regularization is not applicable here, since the gluon propagator behaves as p 2 ~ d and the theory 
remains logarithmically divergent in any dimension d. To our knowledge most other invari- 
ant regularizations such as higher derivatives spoil the masslessness of the propagators, which 
makes higher order calculations virtually impossible. Usually in 1/N calculations the analytical 
regularization of [^] is used. However, this breaks gauge invariance. 

We now consider how we can reconcile gauge invariance with the calculational benefit of 
using massless propagators. First, we break gauge invariance of ( |TTj ) from the beginning by 
introducing a new coupling, A, for the ghost-gluon vertex in (|lj). For A = 1 one recovers the 
original model but the theory remains renormalizable for arbitrary A as well. The only effect 
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will be that the identity (^) will no longer hold. The bare coupling Ao is connected to the 
renormalized one, A, by 

A = Z\\ = Z±Z X Z 2 X Z Z X \ , (14) 

where the Z\ now depend on A. Let us suppose that we used an invariant method, such as 
regularization by higher derivatives [p6| , to regularize this extended theory. Then it immediately 
follows from ( |l3| ) and fll4| ) that the equality A = 1 for the renormalized couplings implies that 
Ao = 1 as well. This leads us to the conclusion that A = 1 is a fixed point, /3a (1) = 0. The 
existence of this fixed point is the key point and it does not depend on the regularization used. 
So using any other regularization can only change the position of the fixed point with in general 
A* = 1 + 0(1/Nf). What is important, though, is that the anomalous dimensions calculated 
at the critical point, 7(A*), are scheme independent and, hence, coincide with the anomalous 
dimensions deduced in the original model (|LT|) . Therefore one can use the regularization which 
is most convenient from the computational point of view. Moreover, since we do not need to 
consider diagrams with external ghost legs, then the only diagrams depending on A are those 
with a ghost loop. As is evident from counting powers of 1/iVy these are themselves 0(1/Nj). 
So at this order it is sufficient to set A = 1. 



In what follows we shall use the A-regularization of [12, 13 j. The propagators of the gluon, 
quark and ghost fields are obtained from (|IT|) as 



xab C / ii \ j fi 1 

Gfxip) = -(^T^ + ^Ja) , D^ V ) = f 2 , D C { P ) = - 2 , (15) 

where we define n to be the combination which arises naturally in the calculations as 

n = N f T F Tr spinor / . (16) 
The value of the amplidude G is derived from (|ll|) as 

G = '" ^w-,) • < 17 > 

where we now set \i = d/2. The new gauge parameter £ is proportional to £ but in what follows 
we shall omit the tilde. The regularization of the theory is carried out by shifting the index of 
the gluon propagator, by taking 

G a AiP) = ^ ^Za (*k + ^Ja) , (18) 

where the factor M 2A is introduced to preserve the canonical dimension of the propagator. 
The divergences appearing in the diagrams as poles in A are removed by the renormalization 
procedure. Throughout we will adopt the minimal subtraction scheme. The technical details 
and subtleties which are inherent to this regularization can be found in the detailed discussion 
of p. ||]. However, we record that up to 0(1/ Nf ) there exists a simple algorithm for the 
calculation of the anomalous dimensions, (27|. To write down the corresponding formula in a 
compact form we introduce a factor u for each the gluon propagator by setting 

G vX -> uG vX , (19) 

so that each diagram with fe-internal gluon lines acquires a factor u k . Then the anomalous 
dimensions of the basic fields $ = {tp, A„, c} can be expressed via the renormalization constants 
Z$ where $o = Z$*&, as follows 

7$ = - 2ud u z£ ] , (20) 

u=l 
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Figure 1: Diagrams contributing to the computation of 7/2. The first graph represents the gluon 
self energy diagrams of Fig. [2|. 

where 

= 1 + z£VA + 4 2) /A 2 + O (-^\ . (21) 

The matrix of anomalous dimensions of the system of composite operators, {O?}, which mix 
under renormalization, is given by 

7i fc = 2ud u z\^ + 5 ifc n fci $7$ , (22) 

U=l 

where the mixing matrix, is defined in the standard way 

Of = Z ik O k , (23) 

so that all Green functions of the operator Of are finite. Again, Z-^ is the coefficient of the 
simple pole in A. The derivation of (p0|) and ( |22"| ) can be found in [27|. Though it should be 
stressed that these formulas are valid only up to 0(1/Nj). We note also the obvious resemblance 



of the formulas (20) and (|22|) with those used in dimensional regularization. 



Before proceeding to the results we underline an advantage of the approach of 27] in 
comparison with the method of the self-consistency equations, (SE), of |l^, [l6| which has been 
more widely used in the 1/iV computations. In the SE method to find the critical exponent 77 
one has to calculate the corresponding renormalized Green functions and then solve the self- 
consistency equations. By contrast to calculate critical exponents using ( |20[) and ( |22|) one only 
needs to know the divergent part of the corresponding diagrams. This is much simpler from a 
computational point of view. As is well known progress in multiloop perturbative calculations 
relies heavily on the possibility of expressing the RG functions through the renormalization 
constants. For calculations of the latter there are various calculational shortcuts. We now 
demonstrate the effectiveness of this approach by calculating the exponent 7/ = 2^(g^) and the 
anomalous dimensions of the [tpifi] operator. 

The diagrams which contribute to Z^, at 0(1 /Nj) are illustrated in Figs, [I] and [2] where in 
the latter the dashed line corresponds to the ghost field. We now take A = 1. Having calculated 
the contributions to Z^ from the diagrams in both figures, we find the following expression for 
77 = 27(3*) at 0(1/ Nj). If we set 

rj = rfr/n + 7/2/n 2 + 0(l/n 3 ) , (24) 
m = C FV0 , (25) 
m = Clrja + C F C AVb , (26) 
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Figure 2: The diagrams contributing to the gluon self-energy at 0(1/Nj). 



where 



then in the Landau gauge 
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The functions <3? and are defined as 

*0i) = ^(2M - 3) + ^(3 - M) - 
$(/x) = /(2m -3) - ^(3-/i) 



^(i) - 



^(/i - 1) , 
i) + v'(i) 



(30) 



where ^(a;) = (In T(x))'. The technical details of the calculations and the values for the individual 



graphs will be given elsewhere [28|. 

Next to determine the quark mass anomalous dimension at 0(1/Nj) one computes the 
anomalous dimension of the associated composite operator [ipip]. At 0(1/Nf) this was derived 



in 1 25] as well as the anomalous dimensions of the basic fields. The diagrams contributing to its 
renormalization constant, Z^, are obtained from those for the quark propagator of Figs. || and 
^ by the insertion of the operator in the fermion lines connected to external vertices. In fact, 
it is sufficient to calculate the diagrams arising from the two 2-loop diagrams in the first line of 
Fig. |l]. The contributions to Zj^ from the other diagrams can be related to the contributions 
of the corresponding diagrams to Z^ by 
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Figure 3: External momenta routing in the quark 2-point function with an operator insertion. 



and 

7m = V + 7^ • ( 32 ) 

This follows if we take the flow of the external momenta p in the propagator diagrams along the 
fermion lines and differentiate the corresponding integrals with respect to p^. Since 



d $ 



Op 2 ) 2 



(33) 



the resulting diagrams will have the same topology as those for the [tpip] operator with the only 
difference being the presence of a new insertion containing 7^. Further, we use the fact that the 
pole part of the diagrams after the subtraction of the divergent subgraphs does not depend on 
the external momenta. So we choose the momenta flow as shown in Fig. ||. The big black dot 
in Fig. |H denotes the insertion of the unit matrix / for the mass operator diagrams, and 7^ for 
the propagator ones. Since the insertion of I or 7^ in the vertex influences only the 7 matrix 
structure of the diagram and does not touch the momentum integral, it can be easily checked 
that the identity ( |3l| ) holds irrespective of the explicit structure of the coloured block on the 
Fig. |3[ Eventually, using ( p2| ) we find for 7 m (g*) 



lm{g*) = Im.i/n + 7 m , 2 /n 2 + 0(l/n 3 ) 



where 



7m, 1 
7m,2 



and 7 a and jb are given by 
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(34) 

(35) 
(36) 

(37) 
(38) 



where the former was determined in JET]]. Having derived a ci-dimensional expression for the 
quark mass dimension, we make several remarks on its properties. First, if one sets ijl = 2 — e 
and expands in powers of e, the coefficients to and including 0(e 4 ) agree with those determined 
from the explicit MS four loop perturbative mass dimension of 0, ||, ^, 1C, 11]. Given that we 



have only evaluated three loop diagrams agreement at this order represents a non-trivial check 
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on our analysis. In addition we have also computed the exponent in an arbitary covariant gauge 
and checked that in the final sum the gauge parameter cancels. These details are presented 
in p£|| . Therefore, we are confident that our value for ^ m .2 is correct. This allows us now to 
produce new information on the mass anomalous dimension at 5-loops in MS. First, we write 
the 0(1/Nj) form of the mass anomalous dimension as 



7m (a) 



- ^C F a + 

oo 

+ ( m r0T, 

r=3 



T F Nf 



-C F 
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C F a r 




(39) 



where the order symbol means that we are ignoring contributions from and six loop terms. 
To extract the coefficients 77150 and 7/451 requires a* at 0(e 5 ) and 0(1/Nj). Although the full 5- 
loop MS QCD /3-function is not yet available, the critical coupling which we require at 0(1/Nj?) 
has already been computed in [ff8|. Hence, we find that 
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where C( n ) is the Riemann zeta function. 



Whilst this gives an indication of the exact form of these coefficients at 5-loops, we have 
another motivation for computing it. Recently, there has been activity in trying to estimate 
the higher order coefficients of perturbative functions in four dimensional field theories from 
knowledge of the known lowest orders and the asymptotic behaviour at high orders, |2^, |3C( |. 
This technique known generally as the asymptotic Pade approximant method, (APAP), has had 
varying degree of success. For instance, the 4-loop /3-function of M = 1 supersymmetric QCD 
had been determined analytically in |3l| up to one unknown parameter. By applying APAP 
methods and its refinement, WAPAP, this parameter was determined numerically. Subsequent 
explicit calculations in [32] produced this parameter analytically and it was in good agreement 
with the APAP estimate for it. On the other hand the original application of the APAP method 
to the 4-loop MS QCD /3-function, |3fj], did not yield as accurate a prediction. In this case the 
coefficients of the polynomial in Nt which appears at 4-loops were estimated with the leading 
coefficient fixed from the known result of the 1/Nt expansion, [Jill]. One of the reasons for a less 
accurate prediction was the appearance of new colour group Casimirs at four loops which were 
absent in lower orders and which therefore had not been built into the estimating procedure. 
Taking account of these issues the method was refined in [33 when the full 4-loop MS result 
became available and new estimates were provided instead for the 5-loop /3-function and quark 
mass dimension. In particular the coefficients of the A^-polynomials were estimated. Whilst 
we believe it will be some time before the full five loop MS renormalization of QCD will be 
performed, if the APAP method for higher order estimation is to serve any useful purpose in the 
interim it is important to determine its reliability. The determination of (^l]) for QCD therefore 
provides us with that test. In [^] the coefficient of the dominant Nf term in the 5-loop 
polynomial was fixed by the leading order large Nf coefficient, m^Q. The results were quoted for 
various values of A^ for two cases. One was where the effect of the new 4-loop colour Casimirs, 
Q4, was accounted for and the other was the case where their presence was ignored. We have 



reproduced the values of these coefficients for the choices of A^ given in [33| together with the 
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numerical values of (|4l|). In order to facilitate comparison we note the relation between our 
notation and that of |33| is 

E 4 = - \T F C F m b0 , £> 4 = - \1%C F m hX . (42) 
From the table it would appear that the WAPAP estimates are not in agreement with fl4l]). 



N c 


Z?4 (w Qa) 


L>4 (w/o Q4) 


£)exact 


2 


8.12 x 10" a 


8.88 x 10" 3 


0.0396 


3 


0.037 


0.037 


0.1083 


4 


0.0891 


0.0831 


0.2049 


5 


0.165 


0.148 


0.3292 


20 


4.31 


3.48 


5.5113 



Table 1. Comparison of APAP results for D4 with the numerical value of the exact coefficient. 

However, it is worth putting the estimates in the context of the other coefficients. For N c = 3, 
for instance, the leading coefficient of the Nf polynomial is O(10~ 5 ) whilst the constant term 
has a WAPAP prediction of 530 which is roughly seven orders of magnitude larger. That the 
WAPAP prediction for 777.51 is within an order of magnitude is perhaps remarkable particularly 
given the nature of the estimation method. By fitting to the large order asymptotic behaviour 
of the perturbation series one is essentially ensuring that the constant term of the polynomial in 
A^ is close to the correct value in the WAPAP as this will always be the dominant contribution 
for a range of values of Nf where Nf is relatively small. Therefore, one would not fully expect 
that the first few coefficients to be reliably predicted. Moreover, in light of our exact evaluation 
of 777.51, it would seem worthwhile to refine the WAPAP estimate of the other coefficients by 
taking our value either as an extra normalization or another constraining number to fit to. 

To conclude our article, we emphasise that we have computed the critical exponents corre- 
sponding to the wave function and mass anomalous dimension in d-dimensions at a new order 
in the 1/Nf expansion in QCD. Not only have we provided new information on the perturbative 
structure of these RG functions but we have also demonstrated the viability of the large A^ 
procedure to compute information on the other important quantities such as those which relate 
to operators in deep inelastic scattering. 
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